Approximating the Tutte polynomial of a binary 
matroid and other related combinatorial polynomials^ 

Leslie Ann Goldberg*^, Mark Jerrum'^ 

"■Department of Computer Science, University of Liverpool, Ashton Building, Liverpool L69 

3BX, United Kingdom. 
^School of Mathematical Sciences 
Queen Mary, University of London, Mile End Road, London El 4^S, United Kingdom. 



Abstract 

We consider the problem of approximating certain combinatorial polynomials. 
First, we consider the problem of approximating the Tutte polynomial of a bi- 
nary matroid with parameters q > 2 and 7. (Relative to the classical {x,y) 
parameterisation, q = (x — l){y — 1) and 'y — y — 1.) A graph is a special case of 
a binary matroid, so earlier work by the authors shows inapproximability (sub- 
ject to certain complexity assumptions) for g > 2, apart from the trivial case 
7 = 0. The situation for q = 2 is different. Previous results for graphs imply 
inapproximability in the region —2 < 7 < 0, apart from at two "special points" 
where the polynomial can be computed exactly in polynomial time. For binary 
matroids, we extend this result by showing (i) there is no FPRAS in the region 
7 < —2 unless NP — RP, and (ii) in the region 7 > 0, the approximation prob- 
lem is hard for the complexity class ^RHIIi under approximation-preserving 
(AP) reducibility. The latter result indicates a gap in approximation complexity 
aX q = 2: whereas an FPRAS is known in the graphical case, there can be none 
in the binary matroid case, unless there is an FPRAS for all of ^RHIIi. The 
result also implies that it is computationally difficult to approximate the weight 
enumerator of a binary linear code, apart from at the special weights at which 
the problem is exactly solvable in polynomial time. As a consequence, we show 
that approximating the cycle index polynomial of a permutation group is hard 
for ^j^RHHi under AP-reducibility, partially resolving a question that we first 
posed in 1992. 

Keywords: Tutte polynomial of a binary matroid, weight enumerator of a 
binary linear code, cycle index polynomial. 
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1. Introduction 



The multivariate Tutte polynomial (in q ^ and 7) of a matroid AA with 
ground set E and rank function tm is defined as follows (see [13, (1-3)]) 

^Tuttc(X;<z,7) = E n >' (1) 

where 7 = {7e}eG-E- The tilde in ^xutte is for consistency with [l^; in general 
we follow the notation used there, so we can conveniently access various useful 
identities. 

An important class of matroids are the graphic matroids, i.e., those that arise 
as the cycle matroid M{G) of some graph G (see Section [TT^ for details). The 
Tutte polynomial of graphic matroids in particular has received much attention. 
For convenience, we offen speak about the Tutte polynomial of the graph G 
rather than the Tutte polynomial of its cycle matroid M.{G). Indeed, this 
polynomial was first defined for graphs, and only later generalised to matroids. 

The Tutte polynomial encodes a large quantity of combinatorial information 



about the matroid [14l |19|. l21l[and the complexity of computing the polynomial 
has been much studied [l|, la 0, 13, 13, HO] • One important motivation for this 
study of the Tutte polynomial is that it includes as a special case (when g is a 
positive integer) the problem of computing the partition function of the Potts 
model in statistical physics. To be more precise about the computational task, 
parameters q and 7 are fixed, and the problem is to compute 2'xuttc(-^; 9, 7) 
for an input matroid A4, where 7 is the constant function with 7e = 7 for 
every ground set element e. In order to introduce the topic, we assume in this 
introduction that the parameters q and 7 are rational, though we shall see below 
that this can be generalised. 

Jaeger, Vertigan and Welsh [3] investigated the complexity of exactly com- 
puting ZTutte(A^; Q, 7) given an input matroid A^. They used a different pa- 
rameterisation of the polynomial, but the problem that they studied is identical 
to the one that we describe here. For the record, they Q (2.2)] define 

T{M;x,y)^ (a;-l)''^-^(^)-''^^^^(2/- 1)1^1-''^'^). (2) 

ACE 

Now substitute ^ — y — 1 and q ~ {x — l){y — 1). Let 7 be the constant function 
with 7e = 7 for every e <E E. Assuming that q (hence 7 7^ 0), Equation ([1]) 
gives 

TiM;x,y) = {x-IY-^'^^^Zt^.UM; {x^l){y-l),j) = (g/7)'--^(^)^Tuttc(A^; 9, 7)- 

(3) 

Unfortunately, any two-parameter version of the Tutte polynomial will omit 
some points. On the one hand, setting y = 1 in ([2]) forces q = {x — l){y — 1) 
to be but setting 7e = 7 = y — 1 = in ([T]) does not force q to be 0. On the 
other hand, the single point g = 0, 7 = in ^ corresponds to an entire line in 
the (x, ?/)-coordinate system, where y = 1 but x can have any value. For this 
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reason, it is sometimes convenient [19|, §2.3] to treat the q ^ case as a limit 
case. We will not need to do this here. 

Jaeger et al. showed that, even when the input is restricted to be a graphic 
matroid, exact evaluation is #P-hard, apart from when q = I and at four 
"special points" . The first three of these are {q, 7) = (4, —2), (2, —2) and (2, —1). 
The fourth is the point {x,y) = (1,1) for which Equation ^ is invalid due to 
division by — evaluation at this point corresponding to counting spanning 
trees in the input graph. As already noted, the line 7 = is also easy in our 
parameterisation. Thus, exactly evaluating the Tutte polynomial of a graph is 
#P-hard, apart from when q ~ 1 and at these special points. Jaeger et al. also 
considered the case in which the parameters are complex numbers, where there 
are additional special points, but we do not consider this case here. As they 
noted, exact evaluation can be done in polynomial time for q — 1 and it can 
be done in polynomial time at some of the special points for large classes of 
matroids. For binary matroids, which are a focus of this paper, the first three 
special points at least are polynomial-time computable. A definition of binary 
matroid is given in Section [131 

Our earlier work [lo| investigates the complexity of approximately com- 
puting ZTuttc{-M]q,f) when M is restricted to be graphic. We are interested 
in determining for which points {q, 7) there is a fully polynomial randomised ap- 
proximation scheme (FPRAS) for the polynomial. An FPRAS is a polynomial- 
time randomised approximation algorithm achieving arbitrarily small relative 
error. A precise definition is provided in Section [T31 We survey the main results 
now, partly because we build on them in this article, and partly to highlight the 
differences in computational complexity between the graphic and binary cases. 

For q > 2 we gave inapproximability results both for 7 < and for 7 > 0. 
As already noted, the case 7 = is trivial. In the "antiferromagnetic" case 
7 < 0, we showed [1] that, apart from at the special point (q — A,"f — —2), 
there is no FPRAS for approximately evaluating the Tutte polynomial of a 
graph unless NP = RP. In the "ferromagnetic" case 7 > 0, we showed [l^l that 
the approximation problem is hard for the logically-defined complexity class 
^RHIIi under approximation-preserving "AP-reductions" . 

The complexity class ^RHIIi of counting problems was introduced by Dyer, 
Goldberg, Greenhill and Jerrum Q as a means to classify a wide class of approx- 
imate counting problems that were previously of indeterminate computational 
complexity. The problems in ^RHIIi are those that can be expressed in terms 
of counting the number of models of a logical formula from a certain syntacti- 
cally restricted class which is also known as "restricted Krom SNP" . ^RHIIi 
has a completeness class with respect to AP-reductions which includes a wide 
range of natural counting problems — see Section [T3] for some examples. Either 
all of these problems admit an FPRAS, or none do. No FPRAS is known for 
any of them at the time of writing, despite much effort having been expended 
on finding one. We conjecture that none exists. Proving counting problems to 
be hard for ^RHIIi with respect to AP-reductions is similar to working with 
the Unique Games Conjecture in the area of approximation algorithms for op- 
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timisation problems, or employing the class PPAD in analysing the complexity 
of Nash equilibria. Since a graphical matroid is a binary matroid, both of the 
hardness results for q > 2 mentioned earlier (for 7 < and for 7 > 0) extend to 
the class of binary matroids. 

The paper [§| also includes hardness results for g < 2 which extend to the 
binary matroid case. For example, there is no FPRAS unless NP = RP if either 
7 or (7/7 is less than —2. The interested reader is referred to [§]. 

The situation is different for q = 2. In this case, we showed Q that in 
the region — 2 < 7 < (apart from at special points) there is no FPRAS for 
approximately evaluating the Tutte polynomial of a graph unless NP — RP. 
However, the most that is known for 7 < —2 (see Q) is that the problem is as 
difficult as approximately counting perfect matchings in a graph, a well-known 
open problem. For 7 > 0, the Tutte polynomial of a graph can be approximated 
efficiently — Jerrum and Sinclair have given an FPRAS [lij . In this paper, 
we show that the problem of approximating the Tutte polynomial of a binary 
matroid is apparently harder. In particular, we show in Theorem [T] that there 
is no FPRAS in the region 7 < — 2 unless NP = RP and that the problem is 
hard for ^RHIIi with respect to AP-reductions for 7 > 0. 

The results in Theorem [T] have interesting consequences for the problem of 
approximating related polynomials. It is well-known that the Tutte polynomial 
of a binary matroid contains as a special case the weight enumerator of a binary 
linear code, which will be defined in Section [31 Hence, we immediately get a 
complexity classification (Corollary [7]) for this problem. This, in turn, allows 
us to make progress on a long-standing open problem about the complexity 
of approximating the cycle index polynomial of a permutation group (see Sec- 
tion We had previously shown Q that there is no FPRAS for this problem, 
unless NP = RP, if the parameter, x, is a non-integer. Using our result about 
the weight enumerator of a binary linear code, we show that the cycle index 
polynomial is as difficult to approximate as #RHHi for every parameter value 
X > 1 (Corollary ini). As we will explain in Section 21 it is at the integer points 
that the cycle index polynomial has combinatorial meaning. 

1.1. Matroid preliminaries 

A matroid M is a combinatorial structure defined by a set E (the "ground 
set" ) together with a "rank function" r^vi : E ^ N. The rank function satisfies 
the following conditions for aU subsets A,B C E: (i) < rMi^) < |^|, (ii) A C 
B implies rj^{A) < rj^iB) (monotonicity), and (iii) rM{A-UB) + rM{Ar\B) < 
r_M{A) + r_\4{B) (submodularity) . 

A subset ACE satisfying rj^iA) ~ \A\ is said to be independent. Every 
other subset A <Z E is said to be dependent. A maximal (with respect to 
inclusion) independent set is a basis, and a minimal dependent set is a circuit. 
A circuit with one element is a loop. 

Suppose that is a matroid with ground set E. Then At is associated 
with a dual matroid A4* with the same ground set E but rank function r_A4* 
given by tm' (^) = 1^1 + {E — A) — rM {E) ■ A cocircuit in is a set that is 
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a circuit in M*; equivalently, a cocircuit is a minimal set that intersects every 
basis. A cocircuit with one element is a coloop. 

We will use the matroid operations contraction and deletion. Suppose e S 
£^ is a member of the ground set of matroid M . The contraction M./e of e 
from Ai is the matroid on ground set E — {e} with rank function given by 
rM/e{A) ^ rM{AU {e}) - rM{{e}), for all A C E - {e}. The deletion M\e 
of {e} from A4 is the matroid on ground set E — {e} with rank function given 
by rM\f.(A) = rj^{A), for all ^ C — {e}. We refer the reader to Oxley's 
book [l^ for a thorough exposition of matroid theory. 

1.2. The Tutte polynomial of a binary matroid 

Let M be a matrix over a field F with row set V and column set E. M 
"represents" a matroid with ground set E. The rank r^(^) of a set of 
columns A in the matroid is defined to be the rank of the submatrix consisting 
of those columns. It is easy to see (see [l^) that a rank function defined in this 
way satisfies the three conditions (i)-(iii) for a matroid rank function presented 
in the previous subsection. Therefore, a set A E is dependent in the matroid 
if and only if the columns in A are linearly dependent as vectors. A matroid is 
said to be representable over the field F if it can be represented in this way. It 
is said to be binary if it is representable over GF(2). 

The cycle matroid of an undirected graph G — {V,E) is the matroid M{G) 
represented by the vertex-edge incidence matrix M of G. In this case, rM{G) (^) — 
\V\ — k{V, A) where k{V, A) is the number of connected components of the graph 
{V, A). We simphfy notation by writing ZTutto(G; q, 7) in place of ZTutto(-M(G); q, 7). 
Since the Tutte polynomial of a binary matroid generalises the Tutte polyno- 
mial of a graph, any hardness result for the latter immediately translates to 
the former. In this context, it should be noted that there is a slight mismatch 
between the definition of the Tutte polynomial given here, and the one used 
in the papers we cite, e.g., 13, 3|. There, the Tutte polynomial of a graph is 



defined using the "random cluster" formulation: 

^Tutte(G;g,7) = I^^""'^^ n > = 9"'' ^Tuttc(G;<z,7). (4) 

ACE eeA 

(Note the absence of a tilde!) But since the two formulations differ only by 
an easily-computable factor g'^', all complexity results, whether about approx- 
imate or exact computation, translate directly. 

For fixed real numbers q and 7 we define the following computational prob- 
lem, which is parameterised by q and 7. 

Problem BinaryMatroidTutte((7, 7). 

Instance A matrix M over GF(2) with rows V and columns E representing a 
binary matroid M. 

Output ZTuttc(A^; 9, 7), where 7 is the constant function with 7e = 7 for 
every e ^ E. 
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1.3. Standard definitions: approximation schemes and approximation-preserving 
reductions 

We are interested in the complexity of approximately solving the problem 
BiNARYMATROiDTuTTE(g, 7). We Start with the relevant definitions. The 
reader who is already familiar with the complexity of approximate counting can 
skip this section. We use the presentation from [10|. 

A randomised approximation scheme is an algorithm for approximately com- 
puting the value of a function / : S* ^ R. The approximation scheme has a 
parameter e > which specifies the error tolerance. A randomised approxi- 
mation scheme for / is a randomised algorithm that takes as input an instance 
a: e E* (e.g., for the problem BiNARYMATROiDTuTTE(g, 7), the input would be 
a matrix M over GF(2) representing a binary matroid JV() and a rational error 
tolerance e > 0, and outputs a rational number z (a random variable depending 
on the "coin tosses" made by the algorithm) such that, for every instance a;, 
Pr [e~^/(a:) < z < e^ J{x)\ > |. The randomised approximation scheme is 
said to be a fully polynomial randomised approximation scheme, or FPRAS, if it 
runs in time bounded by a polynomial in |a;| and £~^. As in [1Q|, we say that a 
real number z is efficiently approximable if there is an FPRAS for the constant 
function f{x) — z. 

Our main tool for understanding the relative difficulty of approximation 
counting problems is approximation-preserving reductions. We use Dyer, Gold- 
berg, Greenhill and Jerruni's notion of approximation-preserving reduction 
Suppose that / and g are functions from S* to M. An "approximation-preserving 
reduction" from f to g gives a way to turn an FPRAS for g into an FPRAS 
for /. Here is the definition. An approximation-preserving reduction from / 
to 5 is a randomised algorithm A for computing / using an oracle for g. The 
algorithm A takes as input a pair (x,e) S E* x (0,1), and satisfies the fol- 
lowing three conditions: (i) every oracle call made by A is of the form (w,S), 
where w G E* is an instance of g, and < (5 < 1 is an error bound satisfying 
i5~^ < poly(|x|, e~^); (ii) the algorithm A meets the specification for being a 
randomised approximation scheme for / (as described above) whenever the ora- 
cle meets the specification for being a randomised approximation scheme for g; 
and (iii) the run-time of A is polynomial in |a;| and e^^. 

If an approximation-preserving reduction from f to g exists we write / <ap 
g, and say that / is AP-reducible to g. Note that if / <ap 9 and g has an 
FPRAS then / has an FPRAS. (The definition of AP-reduction was chosen 
to make this true). If / <ap 9 and g <ap / then we say that / and g are 
AP-interreducible, and write / =ap 9- A word of warning about terminology: 
Subsequent to @], the notation <ap has been used to denote a different type 
of approximation-preserving reduction which applies to optimisation problems. 
We will not study optimisation problems in this paper, so hopefully this will 
not cause confusion. 

Dyer et al. Q studied counting problems in #P and identified three classes 
of counting problems that are interreducible under approximation-preserving 
reductions. The first class, containing the problems that admit an FPRAS, are 
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trivially AP-interreducible since all the work can be embedded into the reduc- 
tion (which declines to use the oracle). The second class is the set of problems 
that are AP-interreducible with t^Sat, the problem of counting satisfying as- 
signments to a Boolean formula in CNF. Zuckerman has shown that #Sat 
cannot have an FPRAS unless RP = NP. The same is obviously true of any 
problem to which #Sat is AP-reducible. 

The third class appears to be of intermediate complexity. It contains all 
of the counting problems expressible in a certain logically-defined complexity 
class, #RHni. Typical complete problems include counting the downsets in a 
partially ordered set Q , computing the partition function of the ferromagnetic 
Ising model with local external magnetic fields 0] , and counting the independent 
sets in a bipartite graph, which is defined as follows. 

Problem #BIS. 

Instance A bipartite graph B. 

Output The number of independent sets in B. 

In it was shown that #BIS is complete for the logically-defined com- 
plexity class ^RHIIi with respect to approximation-preserving reductions. We 
conjecture that there is no FPRAS for #BIS. 

2. Approximating the Tutte polynomial of a binary matroid 

This section provides the proof of the following Theorem. 
Theorem 1. Suppose that q > 2 and 7 are efficiently approximahle. 

1. If^ — 0, orif{q,j) is one of the special points {A, —2) , (2,-2) or (2, — 1), 
then BlNARYMATROlDTuTTE(g, 7) can be solved exactly in polynomial 
time. 

2. Otherwise, ifj < then there is no FPRAS for BiNARYMATROiDTuTTE(g, 7) 
unless NP = RP. 

3. //7 > then #BIS <ap BinaryMatroidTutte((7, 7). 

Some of the parts of the theorem follow from our earlier work in 8] and [13] • 
The main new result is Item ([3]). Its proof follows from 

1. the AP-reduction from #BIS to UniformHyperTutte(2, 1) from our 
paper [lo| . 

2. an AP-reduction from UniformHyperTutte(2, 1) to the problem of 
computing the Tutte polynomial of a binary matroid, where the values 7e 
depend on the input (Lemma [3]), and 

3. implementation of these values 7e using series-parallel extensions on binary 
matroids (Lemma [S]). 

The details are in the following sections. 
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2.1. The Tutte polynomial of a uniform hypergraph 

We have seen one possible generalisation of the Tutte polynomial of a graph, 
namely to binary matroids. Another natural generalisation is to hypergraphs. 
The two generalisation are different, but the relationship between them is inter- 
esting, and and will be exploited in one of our reductions. 

It is typical to define a hypergraph as a pair (V, £) in which V is a set 
of vertices, and f is a set of non-empty subsets of V, called hyperedges. For 
our work on the Tutte polynomial, it will be more convenient to extend this 
definition. Thus, we will use the term "hypergraph" to refer to a pair (V, £) in 
which V is a set of vertices, and f is a multiset of non-empty subsets of V, called 
hyperedges. The reason that the collection £ of hyperedges is a multiset, rather 
than a set, is that it is useful for the Tutte polynomial to allow "parallel" edges 
so that certain operations, such as parallel extensions, which we shall define 
below, can be freely applied fl9'|. A hypergraph is uniform if all hyperedges 
have the same cardinality. 

Let H — (V, £) be a hypergraph. The multivariate Tutte polynomial of H 
was studied (under a different name) by Grimmett jl2| . A definition can be 
found, for example, in [lo| . In this paper we will use the Potts model version. 
Suppose that q is a positive integer and that 7 — {"/f}fe£. Let 

Zpotts(i?;g,7)= E n + (5) 

a:V^[q] /Gf 

where [q] = {0, . . . , g — 1} is a set of q spins or colours, and S{S) is 1 if its argu- 
ment is a singleton and otherwise. Identity ([5]) extends the Tutte polynomial 
from graphs to hypergraphs, but only for positive integer q. It is possible to 
provide a formulation for general q along the lines of ^ , but this is not needed 
in what follows. 

We consider the following computational problem. 

Problem UNiFORMHYPERTuTTE(q, 7). 

Instance A uniform hypergraph H = {V,£). 

Output Zpotts{H', q, 7), where 7 is the constant function with 7^- = 7 for every 

The will use the following lemma, which is an easy consequence of the results 
of 0. 

Lemma 2. #BIS <ap UniformHyperTutte(2, 1). 

Proof. This follows from Observation 2, Lemma 14 and Lemma 15 of fl^. We 
note that [l^l stated a more general definition of "hypergraph" in which hy- 
peredges were taken to be multisets, rather than sets. Nevertheless, the con- 
struction in Lemma 15 actually produces a hypergraph that conforms to the 
definition that we use here. □ 
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2.2. A Potts model characterisation 

Just as we used the Potts model version of the multivariate Tutte polynomial 
of a hypergraph, it will be helpful to have a representation of ^Tuttc(-^; Qjjx) 
in terms of the (multivariate) partition function of the Potts model. See [19l 
Theorem 3.1]. Let be a matroid represented by a matrix M over GF((7) with 
rows V and columns E. For every column e d E, let Ag be the linear equation 
J2i£V -^i-e'^i''') ~ where the arithmetic is in GF((7). The Potts partition 
function of A4 is defined as follows: 

^Potts(A^;g,7)= E lli^ + leSei^)), (6) 

a:V^lq] eEE 

where 

J. , . 1 1, if CT satisfies Ag, 

Oe((J) = < 

I 0, otherwise. 

Also, let ZiBing{M; 7) be a synonym for Zpotts(A^; 2, 7). The Ising model is the 
special case (7 = 2 of the Potts model. Note that 

^isi„g(X;7)= E E n>'^-('^) 

a:V^{0,l} ACE eeA 
ACE a:V~¥{0.1} eCA 

where 7^ = ne6A7e- The number of configurations cr : y — > {0, 1} for which 
HegA = 1 is the number of solutions to the system of linear equations 

Aa = {Ag I e G A}, which we denote #Ayi. Thus, 

^isi„g(M;7)= E^-4#Aa (7) 

ACE 

= E 7a2I^I-'^^(^) 

ACE 

-2l^lZTuttc(A^;2,7)- 

It is interesting to compare definitions ([5]) and ^ in the case g = 2 to see 
how they both arise as natural generalisations of the Ising partition function of 
a graph. In the classical Ising model on a graph, each edge (u, v) contributes a 
factor depending on whether a{u) and a{v) are equal. If we think of this condi- 
tion as asserting that the edge (u, v) is monochromatic, then the extension ^ 
to hypergraphs is immediate. On the other hand, we can equally think of the 
same condition as asserting a{u) + a{v) = (mod 2), which leads us naturally 
to definition ^ for binary matroids. 

2.3. Reduction from UniformHyperTutte(2, 7) 

Consider the following computational problem, which is similar to Binary 
MatroidTutte((7, 7) except that the weight 7 is part of the input. 
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Problem VARBlNARYMATROIDTUTTE(g') . 

Instance A matrix M over GF(2) with rows V and columns E representing a 
binary matroid M.. A positive integer N , given in unary. 

Output Ztuiio(A^; 9,7), where 7 is the constant function with 7g = 2^/-^^ — 1 
for every e E. 

Lemma 3. UniformHyperTutte(2, 1) <ap VarBinaryMatroidTutte(2). 

Proof. Let H = (V, E) be a f-uniform liypcrgraph, an instance of UniformHyper 
Tutte(2, 1). Without loss of generality, assume t > 2 since the result is imme- 
diate for t = 2 (since a 2-uniform hypergraph is a graph, and a graphic matroid 
is binary). Let n = |V| and m = \£\ and assume that these are sufficiently large. 
Let 1 be the constant function which maps every / e f to 1. By definition, 

Zp,us{H;2,l) = Yl 2-°'^°(-), 

where mono((T) denotes the number of hyperedges f G S that are monochromatic 

in configuration a. 

Let £ be the desired accuracy of the AP-reduction and let 6 = e/(mln2). 
Let N be any positive integer satisfying 

6m^(n + ln(16m)) 

We will construct a n x Nm matrix M so that N and M constitute an input 
to Var,BinaryMatroidTutte(2). The rows of M correspond to the elements 
of V. The matrix contains N columns, fi,...,/^, for each hyperedge / G 
£. When we construct the matrix, we choose each of these columns to be 
the indicator vector for an even-sized subset of /, chosen independently and 
uniformly at random. 

Given our construction, it is easy to see that a configuration a which is 
monochromatic on / will satisfy the equations corresponding to all N columns 
/i) • • • ) /jv- For this, it is important that the random subsets of / corresponding 
to these columns have even size since the relevant equation Aj^. is X^^gy Mij.a{i) = 
mod 2. 

Suppose that a configuration a is not monochromatic on a hyperedge / 
and that it assigns i elements of / to spin 1 and k elements of / to spin 
for positive integers £ and k. Note that £ + k = t since the hypergraph is 
t-uniform. The number of even-sized subsets of / is 2*~^ and the number of 
even-sized subsets for which <t restricted to that subset has an even number 
of Is is 2^~^2'^~^ = 2*~^. Thus, the probability that the equation associated 
with a column /, is satisfied is 1/2. So, by a Chernoff bound, the number 
of columns fi with satisfied equations is with high probability in the range 
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[(A^/2)(l - S), {N/2){1 + 5)]. Specifically, the failure probability for this event 
is at most 



2exp(-52jv/6) < 2cxp(-n - InlGm) < 2 x 2"" x ^ 



1 



16m 8m2"- 



171 . Cor. 4.6]. From the miion bomid — ranging over events indexed by the 2" 
choices for a and m choices of column / — we conclude that, with probability at 
least 7/8, the following is true. For every configuration cr, for every hyperedge 
f & £ on which a is not monochromatic, the number of columns in {/i, . . . , /jv} 
with equations satisfied by a is in the range [(A^/2)(l — S), (7V/2)(1 + S)]. 

Let M be the binary matroid represented by M. Let y = 2^/''^ and 7 = y — 1- 
Let 7 be the constant function which maps every element of the ground set 
of to 7 as in the definition of VarBinaryMatroidTutte. To complete 
the verification of our reduction, we need to show that Zpotts(^^; 2, 1) may be 
easily computed given ZTuttc(A4, 2, 7); note that the latter quantity is equal 
to 2^I^IZising(A^; 7) by (O. Now the contribution of a configuration a to the 
quantity 2™Zpotts(i?; 2, 1) is 2'"2'"™°('") . Let be the contribution of a to 
^ising(X;7)- Then 

< yJVmono(o-)y(m-mono(o-))(Af/2)(l+(5) 

Also, we get a similar lower bound. 

^cr > y^'n™°('^)y(™-™°n°('^))(^/2)(l-'5) 

^ ^mA^(l — (5)/2^(A^/2)mono(CT) 2^^^<^2^2^°^°^^^ 

The reduction has the desired accuracy for an AP-reduction, since 2™* — e^. 
Assume that the oracle caU to VarBinaryMatroidTutte(2) is powered to 
have failure probability |. Then the overall failure probability is bounded by 
J, being the sum of | from the randomised nature of the reduction itself, and 
I from the single oracle call. □ 



2.4- Series-parallel extensions of binary matroids 

The standard method for reducing the problem of evaluating the Tutte poly- 
nomial with some weight 7' to a Tutte-polynomial evaluation problem with a 
different weight 7 is to "implement" the edge weight 7' using series and parallel 
extensions of weight-7 edges. See [l^. Section 2.3] and [13, Section 10]. Series 
and parallel extensions of matroids are generalisations of the stretchings and 
thickenings of graphs used by Jaeger et al. [13. Sokal [H Section 2.3] has given 
the details, both for graphs and for general matroids. It is fairly easy to show 
that these extensions can be done within the class of binary matroids. We do 
this here. 
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Lemma 4 (parallel extension for binary matroids) . Let M. be a binary matroid 
represented by the matrix M with rows V and columns E. Let 7 = {TejesB- 
Let c be any column in E. Suppose 71 > and 72 > satisfy 

l+7c- (l + 7i)(l+72). (8) 

Let e! be a new column. Define 7' = {"f'e}eeEu{e'} follows. Let 7^, = 72 and 
7^ = 71 . For every other column e G E, let 7^ = 7e • There is a binary matroid 
A4' represented by a matrix M' with \V\ rows and columns E U {e'} for which 

ZTuttc{M;q,j) = ZTuttc(A^';g,7')- 

Proof. Let M' be the matrix constructed from M by making column e' a copy of 
column c. Let M-' be the matroid represented by M'. Note that, for any ACE, 
rM{A) = rM'{A). Thus, from the definitions in Section [TTl M = M'\e'. The 
result now follows from [l^, (4.22)], provided that we can show that either (i) 
c and e' form a two-element circuit of Ai', or (ii) c and e' are both loops of 
Ad' . (These are the two side-conditions for the application of [l^ (4-22)].) 
Now, If c is not the all-zero vector, then c and e' do form a two-element circuit 
(minimal dependent set) of M' since rM'{{c, e'}) = rM'{{c}) = rM'iW}) = 1- 
Otherwise, rM'{{c}) ~ rM'iW}) — Oj so both c and e' are loops of M' . □ 

Lemma 5 (series extension for binary matroids). Let M be a binary matroid 
represented by the matrix M with rows V and columns E. Let 7 = {le}eeE- 
Let c be any column in E. Assume q ^ 0, and suppose 71 > and 72 > satisfy 

Let e' be a new column. Define 7' = {l'e}eGEu{e'} o,s follows. Let 7^/ = 72 and 
7^ = 71 . For every other column e ^ E, let 7^ = 7e • There is a binary matroid 
A4' represented by a matrix M' with \V\ + 1 rows and columns E U {e'} for 

which (1 +7i/g + 72/g)ZTutto(A^;g,7) = ^Tutto(A^'; g, 7')- 

Proof. Let M' be the matrix constructed from M by adding a new column e' 
and a new row r' . The new row has ones in column c and column e' only. There 
are no other ones in column e'. Let M.' be the matroid represented by M' . We 
first show that, for any A C E, 

rM{A)^rM'{AUe')~l. (10) 

Equation (|10p can be verified by checking three cases 

• A C {c}, 

• c e A and rM{A^} — fMi^ ~ {c}), and 

• c€ A and rM{A) =rM{A- {c}) + 1. 
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(The conditions on the rank function of a matroid guarantee that these are the 
only three cases.) Now note that rM'iW}) = 1, so (ITOl) imphes that ^^(A) = 
rM'{A U je'}) — rM'{{e'}). We conclude (from the definitions in Section [TTT|) 
that AA = M! I e' . The result now follows from [H (4.28)], provided that we 
can show either (i) c and e' form a cocircuit oi M! or (ii) c and e' are both 
coloops of M.' . (These are the two side-conditions for the apphcation of [l^ 
(4.28)].) 

Suppose first that c is a coloop of M. (i.e., an element that is present in 
every basis oi M). Consider any independent set A C — {c} of M' . Since 
c is linearly independent of the columns in A in M, both c and e' are linearly 
independent of the columns in A in M' . Thus, both c and e' are coloops of M! 
and we have (ii). 

Finally, suppose that c is not a coloop of M.. Our goal is to prove (i). First, 
Any independent set in M including neither c nor e' can be extended to a larger 
independent set by adding either one of c or e', so every basis of M.' intersects 
{c, e'}. However, since c is not a coloop of A^, there is a basis i? of that does 
not include c. Then B U {c} and B U {e'} are both bases of M' and so {c, e'} 
is a minimal set that intersects every basis oi M! . Thus, we have (i). □ 

2.5. Implementing variable weights 

Lemma 6. Suppose that 7 > is efficiently approximable. Then VarBinary 
MatroidTutte(2) <ap BinaryMatroidTutte(2, 7). 

Proof. Let C-^ be a sufficiently large function of the parameter 7. The exact 
computation of C-y is from 10|. This will be explained below. 



Let M and N be an instance of VarBinaryMatroidTutte(2). Let M be 
the matroid represented by M. Suppose that M has n rows and m columns and 
that the product Nm is sufficiently large with respect to the constant Cy. Let 
7' = 2^/^ — 1. Let 7' be the constant function which maps every ground set 
element of to 7'. 

The proof is based on the proof of [10|, Lemma 17]. Let e be the desired 
accuracy in the approximation-preserving reduction. Let 



X 



£2 



Let 7 be a rational in the range e~^7 < 7 < 6^7. Since 7 is efficiently approx- 
imable, the amount of time that it takes to compute 7 is at most a polynomial 
in m, N and e~^. 

The idea of the proof is to show how to use series and parallel extensions of 
weight-7 elements to implement weight 7* satisfying 

e-xy < 7* < 6^7'. (11) 

Let 7* be the constant function which maps every ground set element of A4 to 
7*. The definition of ^Tutto and the fact that x ^ e/(4m) imply that 

e-"/*^Tuttc(A^;2,7') < ZTuttc(X;2,7*) < e^/4^Ti.tte(A^; 2, 7')- 
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Let 7 be the constant function which maps every ground set elemeiit^to 7. 
We can think of our implementations as constructing a binary matroid A4 such 
that ^Tutte(A4; 2,7*) is equal to the product of ^Tutto(-^; 2, 7) and an easily- 
computed function of 7. This easily-computed function arises from the extra 
factor (1 -|- 71/q' -|- 72/9) in Lemma [5] We will ensure that the matroid A4 has 
at most C-yTn^N/e ground set elements. To finish, we note, using the definition 
of 2'Tutto and the definition of x, that 

e-"/^^Tuttc(A?;2,7) < ^Tuttc(>l;2,7) < e^/^ZTuttel^l; 2, 7), 

where 7 is the constant weight function which assigns every element weight 7. 
We finish the approximation of ^Tuttc(-A4; 2,7) by using the oracle to approxi- 
mate 2'Tutto(A^; 2,7) using accuracy parameter 5 = e/2. 
It remains to show how to do the implementation. Take 

^=|(22/^-l) = |7'<7'(l~e->^). 



The proof of 10|, Lemma 17] shows how to use series and parallel extensions 
of weight 7 elements (from Lemmas |4] and [5]) to implement a weight 7* which 
satisfies 7' — tt < 7* < 7'. This ensures that Equation (fTT|) holds. The series and 
parallel extensions in the implementation of 7* introduce at most C^log(7r~^) 
ground set elements, where is some quantity depending on 7 but not on n, m, 
or N. Note that 2^/^ - 1 > 21n(2)/iV so tt"! < 2N/x- Thus, mC^ log(7r-i) < 
C^m?N/e and the matroid Ai has at most C^rr?N je ground set elements, as 
required above. □ 

2. 6. The Proof of Theorem [7] 

Proof of Theorem [71 Item ([T]) comes from 14 1 . Item ([2]) follows quickly from [§| 
together with an application of matroid duality; here are the details. When 
g > 2, item ^ follows from the corresponding hardness result for graphs, since 
binary matroids generalise graphic matroids. The same is true when q — 2 and 
7 > —2. The complementary case q = 2 and 7 < —2 then follows by matroid 
duality, as can be seen by combining the following observations: Binary matroids 
are closed under duality [3, 2.2.9] and the representation of the dual of a binary 
matroid can be constructed efficiently [H, 2.2.8]. Also, if A4* is the dual of M. 
then ZTutte{-M', g, 7) is an easily-computed multiple of .^Tuttc(A^*, 1, l/l) — see 
[H, (4.14)]. Finally for g = 2, if 7 < -2 then -2 < q/j < 0. 

Item ^ follows from [l3, Theorem 1] for q > 2 since the cycle matroid of 
a graph is a binary matroid. For q — 2 it follows from Lemmas [21 Lemma [H 
and [a □ 

3. The weight enumerator of a binary linear code 

Given a generating matrix M over GF(2) with r linearly independent rows 
and c columns, a code word w is any vector in the linear subspace T generated 
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by the rows of M . For any real number A, the weight enumerator of the code 
is given by VFm(A) — X^toex where |jw|| is the number of non-zero entries 

in w. We consider the following computational problem, parameterised by A. 

Problem WE(A). 

Instance A generating matrix M over GF(2). 
Output Wm{X)- 

It is well-known (see below) that the weight enumerator of a binary linear 
code is a special case of the Tutte polynomial of a binary matroid. Thus, 
Theorem [T] has the following corollary. 

Corollary 7. Suppose that A is efficiently approximable. 

1. If X G { — 1, 0, 1} then WE(A) is solvable in polynomial time. 

2. // |A| > 1 then there there is no FPRAS for WE(A) unless NP = RP. 

3. //A e (-1,0) then there there is no FPRAS for WE(A) unless NP = RP. 

4. //A e (0,1) then #BIS <ap WE(A). 

Corollary [7] follows immediately from Theorem [1] and from Lemma [8] be- 
low, which is the q — 2 case of a result of Greene [Til, Corollary 4.5]. See 
also Cameron 0, Theorem 4.1], but note that both authors employ a different 
parameterisation of the Tutte polynomial. We provide a short proof here for 
completeness, since we have already done almost all of the necessary work. 

Lemma 8 (Greene). Let M be a generating matrix over GF(2) with rows V 
and columns E. Let M. be the binary matroid represented by M . Let A be any 
non-zero real number and let ^ = — \. Let 7 be the constant function with 
7g = 7 for every column e of M . Then 

I^m(A) = Al^l2l^lZTuttc(X;2,7)- 
Proof. As in the proof of Lemma [3l we have, from equations (|6]) and ([7]), 

2l^lZTuttc(A^;2,7) = ^isi„g(A^;7)= J2 (l + 7r'^"^= E A-^^^^-^), 

(12) 

where sat(cr) denotes the number of columns e G -B such that 

Y,M^,e<J{i)^0 (mod 2), 
iev 

or, using the notation introduced earlier, (5e((T) = 1. Similarly, 

W^m(A)- J2 n^'"'^^"^= E Al^l-^-*(-), (13) 

<T:V^{aS} eeE a-:V^{QS} 

where, corresponding to cr, the codeword w is the sum of the rows i with a{i) — 1, 
so 1 — (5e(cr) is the bit in position e of the code- word. The result follows by 
comparing and (IT51) . □ 
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4. The cycle index polynomial 

Let r be a group of permutations of {1, . . . ,1/}. Each permutation .g G T 
decomposes the set into a collection of cycles, eye (5) denotes the 

number of cycles in this decomposition. The single- variable cycle index polyno- 
mial of r is defined as follows. 

' ' ger 

When X is a positive integer, Zci(r;a;) counts the number of orbits (or 
equivalence classes) when strings from a size-x alphabet are operated on by 
permutations in F (which permute the character positions in the strings). Two 
strings are in the same equivalence class if there is a permutation in F which 
maps one into the other. For example, when x = 2, v = 3, and F is the 
symmetric group (on 3 elements), the elements of F are the identity permutation 
(which has 3 cycles), the three transpositions (12), (13) and (23) (each of which 
has 2 cycles), and the singleton cycles (1 2 3) and (1 3 2). Thus 

Zci(F; 2) = i (2^ + 3 X 2^ + 2 X 2^) 4. 

Thus, there are four orbits (namely the orbits of the strings 000, 001, Oil, and 
111). For more details, see [1, [13, We consider the following problem, in 
which the parameter x is a positive real number. 

Problem CYCLElNDEX(a;) 

Instance A set of generators for a permutation group F 

Output Zci{T;x). 

We showed @, Theorem 4] that if x is not an integer then there is no FPRAS 
for CYCLElNDEX(a;) unless RP — NP. In fact, it is NP-hard to approximate 
Zci(T]x) within any polynomial factor. However, our technique from The- 
orem 4] does not say anything about the difficulty of the problem in the more 
interesting case when x is an integer. We raised this question in Q but were 
unable to resolve it (until the present paper). Note that it is easy to compute 
Zci(T] 1) exactly in polynomial time. Corollary[7]has the following consequence. 

Corollary 9. Suppose that x > is efficiently approximahle. Then 

1. If X ^ 1 then CycleIndex(x) can be solved exactly in polynomial time. 

2. If X is not an integer then there is no FPRAS for CYCLElNDEx(a;) unless 
NP = RP. 

3. If X > 1 is a positive integer then #BIS <ap CycleIndex(x) . 

Proof. Items ([T]) and ^ are from [gI, Theorem 4] . We now prove item ([3|) . Let 
A = x~^. Note that A e (0,1). We know from Corollary [7] that #BIS <ap 
WE(A). To finish, we show that WE(A) <ap CYCLElNDEx(a;). 
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The reduction is straightforward. Let M be a generating matrix (an instance 
of WE(A) with r hnearly independent rows and c columns). Let T be the 
subspace generated by the rows of M (this is the set of code words). Let Mj_* 
denote row i of M . For each binary string m G {0, 1}'', let w{m) be the vector 
Si=i "m-iMi^^. (where arithmetic is over GF(2)) and let fc(m) be the number of 
I's in this vector. Since the rows of M are linearly independent, each element 
of T is uniquely expressible as w{m) with m G {0, 1}'', so 

Wm{\)= (14) 

me{n,l}'' 

Let V — 2c. Our objective will be to construct a group F of permutations of 
{1, . . . , i^} so that VFm(A) can be computed from Zci(F; x). For i G [r], let Qi be 
the permutation of {1, . . . , i'} defined as follows: For all j G [c], elements 2j — 1 
and 2j are mapped to each other by gi if Mi,j = 1 and each of these is mapped 
to itself by gi if M^.j = 0. Let F be the group of permutations of {1, . . . , 
generated by 51, ... ,5^. 

For each binary string m G {0, 1}'', let g{m) be the permutation g"^^ ■ ■ ■ g™'', 
where g^ denotes the composition of £ copies of the generator gi so g° is the 
identity permutation. Note that, for each j G [c], elements 2j — 1 and 2j are 
swapped by g(m) iff X]I=i "^iMi^j — 1 mod 2. Thus, the number of j for which 
there is a swap is fc(m) so cyc{g{m)) = v — k{m). 

We will show that each permutation g € T can be written as g{m) for ex- 
actly one m G {0, 1}''. First, suppose that g — gii ■ ■ ■ git for some ii, . . . ,ie G 
{1, . . . ,r}^. Since the generators gi,. . . ,gr commute, we can re-order so that 
ii,...,i£ are monotonically non-decreasing. Then, since the generators all 
have order 2, we can cancel factors that are the identity permutation, making 
ii, . . . , if distinct. Thus, g can be written as g(rn) for at least one m G {0, 1}''. 
To see that m is unique, suppose that g{m) = g{rn') for m ^ m' . Then 

ffi"^ ■ • ■ ff™' — ff™^ ' ' ' 9r^ ■ Thus, for each j G [c], the number of permuta- 
tions in I mi = 1} which swap element 2j — 1 with element 2j has the same 
parity as the number of permutations in {(^i | — 1} which swap element 
2j — 1 with element 2j. Thus, 

r r 

TxiiMi^^ — ^ m'^Mi^^ mod 2 

i=l 1=1 

Since the rows of AI are linearly independent, we conclude that m = m' . Thus, 
we have proved that F — {g{m) | m G {0, 1}''} so, from the definition of 

|F|Zci(F;a;) = J2 x"'^'^'''^ ^ x" ^ X''^'^^ ^ x^'WmW, 
where the last equality uses Equation ([T^. □ 
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